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Abstract 

The structure of single atoms in real space is investigated by scanning tunneling microscopy. Very 
high resolution is possible by a dramatic reduction of the tip-sample distance. The instabilities 
which are normally encountered when using small tip-sample distances are avoided by oscillating 
the tip of the scanning tunneling microscope vertically with respect to the sample. The surface 
atoms of Si(lll)-(7x7) with their well-known electronic configuration are used to image individ- 
ual samarium, cobalt, iron and silicon atoms. The resulting images resemble the charge density 
corresponding to 4f, 3d and 3p atomic orbitals. 
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FIG. 1: (Color online) Principle of dynamic STM. The tip is mounted on a quartz cantilever (not 
shown here, see ^(| for details) with a stiffness of 1800 N/m which oscillates at a fixed amplitude 
A ~ 0.3 nm. The unperturbed resonance frequency is /o ~ 20 kHz. Conservative components 
in the tip-sample force are measured by the change in the cantilevers eigenfrequency they are 
causing 18]. Non-conservative components in the tip-sample force are measured by monitoring 
the dissipation energy which has to be supplied to the cantilever to maintain a constant amplitude 

Q. 

Democritos of Abdera (460 BC - 370 BC) predicted the existence of atoms and speculated 
that atoms come in different sizes and shapes yj. Today, we know that the size and shape 
of an atom depends on its chemical species and the way it is bonded to neighboring atoms. 
Atoms consist of the tiny, positively charged nucleus and electrons with negative charge — e 
which orbit the nucleus with a speed of a few percent of the speed of light. Because of the 
uncertainty principle of quantum mechanics, articulated by Werner Heisenberg in 1927, we 
can not observe single electrons within their orbit. However, quantum mechanics specifies 
the probability of finding an electron at position x relative to the nucleus. This probability is 
determined by |?/>(x)| 2 , where ?/>(x) is the wave function of the electron given by Schrodinger's 
equation |2]. The product of — e and |-?/>(x)| 2 is usually interpreted as a charge density, 
because the electrons in an atom move so fast that the forces they exert to other charges 
are essentially equal to the forces caused by a static charge distribution — e|^(x)| 2 . For the 
electrons of an atom, the solutions to Schrodinger's equation are wavefunctions ^^(x) 
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defined by four quantum numbers: principal quantum number n, angular momentum I, 
^-component of angular momentum m and spin a jjj. 

The shape of atoms, as given by their total charge density, can be determined by X-ray- or 
electron-scattering if they are arranged periodically in a crystal. X-ray scattering allows to 
measure the roughly spherical charge density of ions in alkali-halides or the non-spherical 
structure of atoms crystallized in a diamond lattice [5|. Recently, d-orbital holes were 
found in the charge density of copper-oxide compounds by transmission electron microscopy 
measurements . Observing the charge density of single atoms in real space became feasible 
in 1982 with the invention of scanning tunneling microscopy (STM) [2|. In contrast to the 
techniques listed above, STM does not probe the total charge density of an atom, but the 
charge density at the Fermi level p Fermi |3 S|- For negative sample bias voltage, the charge 
density of the energetically highest occupied electron states is imaged, while for positive 
sample bias, the charge density of the lowest unoccupied states is imaged. If the states 
at the Fermi level have an s-type symmetry, the const ant- density surface pFermii*) = c is 
composed of spheres centered at the nuclei of the surface atoms. If p Fermi is built from atomic 
states with higher angular momentum, ppermi can have a more complex shape, in particular 
for large values of p Fermi m close proximity to the atoms. Thus, if higher momentum states 
are present at the Fermi energy, the corresponding STM images are expected to reflect 
their shape, such that the shapes of atomic images are non-spherical and display orbital 
substructures. However, it is important to note that orbitals in the sense of mathematical 
constructs with an amplitude and phase can not be observed, but the density of charges 

nn 

within a certain energy range corresponding to specific orbitals, see [8, 9] and the discussion 
[lol ] and references therein. 

In spite of the large amount of standard STM data gathered, the observation of atomic 
substructures linked to atomic higher-momentum states has not been reported, instead 
STM images of atoms always appear more or less parabolic. Interestingly, tunneling chan- 
nels originating from higher orbital momentum states have been found in tunneling exper- 
iments involving mechanically controllable break junction techniques jljj, where tunneling 
gap widths of the order of bulk next-neighbor distances can be realized. However, imaging 
is not possible with the break junction technique. While 'subatomic' resolution, showing 
two maxima in the i mag e of a single atom linked to sp 3 -type orbitals has been observed by 
force microscopy jl^ . Q], features within single atoms have not yet been observed by STM. 
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After a presentation of the work reported in by a member of our team, G. Doyen from 
the University of Munich insisted that orbital structures should be visible in STM, if these 
structures can be imaged by force microscopy. Also, C.J. Chen (chapter 7 in discusses 
the similarity of the functional dependence of tunneling current and attractive part of the 
tip-sample force. These theoretical arguments and the experimental appearance of higher 
momentum states in AFM- and break-junction-tunneling- experiments inspired us to think 
about possible reasons for the absence of subatomic resolution in standard STM jlHI ]. 

One likely cause is the large tip-sample separation in standard STM (0.5 - 0.8 nm). 
Even if the Fermi level in a sample or tip was occupied with higher momentum states, the 
contour lines of the electron density show at this large distance only one single maximum 
per surface- and tip-atom. A second possible cause is the geometry of the electronic states 
participating in the tunneling process. The valence bands of solids can often be constructed 
from atomic states by Linear Combination of Atomic Orbitals (LCAO). As a consequence, 
spatial features are expected to be present in individual tunneling channels Q| . Because of 
various effects jjj, the energy of atomic (and bulk-) states depends on the angular momentum 
I. Thus, if large tunneling voltages are used, many states can contribute to the tunneling 
current, smearing out spatial features. We concluded that small bias voltages and very small 
tip-sample distances are required for enhanced resolution. 

In standard STM, the tip sample distance decreases when the tunneling impedance, given 
by the bias voltage divided by the setpoint of the tunneling current is decreased. In STM 
operation, the setpoint of the tunneling impedance has to satisfy a compromise. If the 
setpoint of the tunneling current is too small, the corrugation in the image is small and the 
noise level is high. If the current setpoint is too high (and/or bias voltage too low), the 
tip-sample forces become excessive and destruction of tip and sample occurs when scanning 
at high loading forces. However, when the lateral scan is stopped to allow for current-versus- 
distance spectroscopy experiments, small tip-sample distances can be realized without tip 
or surface destruction. A recent theoretical study has shown that single atom tips can 
withstand repulsive forces up to 30 nN without damage. However, pushing the tip laterally 
over a surface at high loading forces damages both tip and sample. A stop-and-go approach 
where the scan is stopped, the tip is lowered until it reaches its high current setpoint, 
retracted and moved laterally to start over again would remedy this situation. We found 
experimentally that a much simpler and faster approach is possible by deliberately vibrating 
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the tip vertically (see Fig. [TJ). Because of the strong distance dependence of the tunneling 
current, the current only flows when the tip is close to its lower turnaround point. The lateral 
motion of the tip during this time when it is very close to the sample is given approximately 
by the scanning speed times less than half of the oscillation period, amounting to only a few 
pm. 

A similar concept already proved successful in tapping force microscopy, where the oscil- 
lation of the tip also reduces lateral forces and decreases tip wear Q]. The tip oscillation 
in dynamic STM reduces or even prevents the wear which normally occurs when scanning 
a probe tip across a surface with repulsive load forces. 

In our experiments, the vibration of the tip is done by mounting it on a quartz force sensor 
which oscillates at frequency fo ~ 20 kHz, as described in detail in Q]. While the oscillation 
could also be fashioned with a standard STM, mounting the tip on a cantilever allows to 
monitor two more important observables: the conservative and the dissipative components 
of the tip-sample force. The forces which occur during tunnelin g ha ve been measured before 
with a different technique, albeit not with atomic resolution 13] . Measuring the conser- 
vative force component yields information about the tip-sample distance. Monitoring the 
dissipative component allows to assess if tip and sample atoms are deflected strongly from 
their equilibrium positions [19|- The STM images are recorded in the topographic mode, 
where the average tunneling current (averaged over a time of approx. 30 oscillation cycles) 
is kept constant by adjusting the ^-position of the tip accordingly. For an exponentially 



ing current I av is smaller 
. In this publication, we 



decaying tunneling current I(z) = Ioexp(—2Kz), the average tunned 
than the peak tunneling current I p k with I av /Ipk ~ (2ttk,A)~ ' 5 
present measurements performed by dynamic STM and report the observation of atomic 
images containing symmetry features linked to p, d and f orbitals of Si, Co and Sm atoms. 
All experiments have been performed in ultrahigh vacuum at a pressure of p ~ 10~ 8 Pa and 
ambient temperature T m 300 K. 

In STM, the role of tip and sample is perfectly interchangeable - the image is a convolution 
of tip and sample states. For zero temperature, the tunneling conductivity depends on the 
matrix elements in a first order calculation as 



£ = (27r) 2 G J2 D Sf,(E F )D Tu (E F )\M^\ 2 (1) 
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with the conductance quantum Go = 2e 2 /h. Here, the densities of states Dsn and Dt u 
of the corresponding band states of sample and tip as well as the matrix elements M^ v are 
assumed to be constant in the energy range eV around the Fermi level Ep. 

In Bardeen's treatment of a tunneling junction [9], the matrix elements are evaluated 
to be 

Mr = I feV 2 ^ - ^W.) dV (2) 

zm e Jn T 

with the electron's mass m e and tip and sample states Xu and ip^, respectively. The 
integration extends over the tip volume Qt, see e.g. chapter 2 in |9j. The tunneling process 
is considered to be spin-independent for the case of a non- magnetic tip and/or sample. The 
modulus of the matrix element does not change if tip and sample states switch roles, giving 
rise to the reciprocity principle: If the electronic states of the tip and the sample under 
observation are interchanged, the image should be the same - see p. 88 in Jj|. Thus a STM 
image can either be interpreted as one image of N sample atoms probed by the front atom 
of the tip or as N images of a front atom probed by iV sample atoms. In this study, we are 
interested in the image of the front atom of the tip. Therefore, the electronic state of the 
sample atoms needs to be known. We have chosen a Si (lll)-(7x7) surface as sample (see 
Fig. El (a), because the electronic states of the surface atoms (adatoms) in Si (lll)-(7x7) 
are precisely characterized, so that the tip states can be inferred from the images. 

Each adatom of the Si (lll)-(7x7) surface exposes one sp 3 -dangling bond oriented per- 
pendicular to the surface. The large distance between adjacent adatoms (at least 0.67 nm) 
allows to clearly separate individual tip atom images. The shape and apex radius of the 
images is used to characterize the electronic states of the tip. The advantage of scanning an 
area containing many sample atoms is that potential errors due to multi-tip effects are easily 
detected. The observation of single defects on the Si surface rules out double- or multi-tip 
effects, and the observation of deep corner holes offers a proof that the tip is pointed and 
sharp. This is not a restriction of our method, because in principle, well defined tips can 
also be created 



1131 ] and unknown sample states could be studied with these tips. As a 
reference, standard STM images of Si (lll)-(7x7) are shown in Fig. [21 Fig. El (a) shows the 
structure of the Si (lll)-(7x7) surface and Figs. El (b), (c) are STM images taken with an 
iridium tip, using a tunneling current of 1 nA and a sample bias of +2.4 V (empty states) 
in Fig. Eb, and a sample bias of -2.4 V (filled states) in Fig. El (c). The apex of the adatom 




FIG. 2: (Color online) (a) Surface unit cell of Si (lll)-(7x7) after the dimer-adatom-stacking fault 
model. The diamond shaped unit cell has characteristic features: deep cornerholes at the corners 
of the diamond (see Fig. 1 for a cross section of the corner hole along the long diagonal of the unit 
cell), 18 dimer atoms, 12 adatoms and 6 restatoms. One triangular halve of the unit cell develops 
a stacking fault. The surface atoms (adatoms) are printed in red. Each adatom is bonded to the 
layer underneath by three covalent bonds. The forth valence electron forms a dangling bond (3sp 3 
state), which sticks out perpendicular to the surface. The rest atoms, printed in full black circles, 
are located more than 100 pm below the adatoms. (b) Standard STM image of the empty states on 
the Si (lll)-(7x7) surface (sample bias 2.4 V, tunneling current 1 nA). (c) Standard STM image 
of the filled states on the Si (lll)-(7x7) surface (sample bias -2.4V, tunneling current 1 nA). In 
the filled state image, the upper half of the diamond-shaped unit cell appears to be slightly higher 
than the lower half, caused by the stacking fault. In images b and c, the image of a single adatom 
is a paraboloid with a radius of ~ 0.5 nm. Ir tips have been used for both images. 

image can always be approximated as a paraboloid. The topographic data around the peak 
is approximated with z(x,y) = z — (x — xp ea k) 2 /2R X — (y — yp ea k) 2 /2R y , thus R is easily 
determined from the topographic data. Application of this formula to the measured profiles 
in Fig. I2t»,c yields apex radii of 0.5 nm. 

Figure 01 shows the difference between imaging in standard STM and dynamic STM. A 
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FIG. 3: (Color online) (a) STM images of a Si (lll)-(7x7) surface taken with a Sm tip. A 
Standard STM, sample bias: 2 V, I t = 100 pA, (b) Standard STM, sample bias: -2 V, I t = 200 
pA, (c) Dynamic STM, sample bias: -1 V, I av = 100 pA, A = 0.8 nm. 

samarium tip was used to image Si in standard STM in Fig. El (a) and Fig. El (b), and dynamic 
STM in Fig. El (c). Because Sm is a relatively soft material, the tip can only withstand a 
small tunneling current of 100 pA at a sample bias of 2V in standard STM. In dynamic 
STM with A = 0.8 nm, the average current can be set to 100 pA (yielding a peak current 
I Pk pa 700 pA) at a bias of -1 V. Thus, dynamic STM allows to get much closer to the surface 
without tip degradation, resulting in improved spatial resolution. The curvature of an STM 
image of a single atom is a characteristic feature. According to a calculation by Chen j^, 
the radii of the atom images is a function of the distance z between tip and sample atoms 
and the type of atomic tip and sample states: 



s — s : R 



(3) 



s - Pz ■ R 



1 + 1/kz 



(4) 
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FIG. 4: (Color online) (a) Apparent radius R of a surface atom in standard STM mode as a 
function of tip and sample states for a tip-sample distance of z = 0.4 nm. Calculated with Eqs. (jSJ- 
EJ) after section 6.2 in L 9j. The tunneling current is assumed to decay with distance d as exp{—2nd) 
with k = 10 nm -1 . a Tip state: s orbital, sample state: s orbital. The apparent radius of the 
adatom image is given by the distance, thus R a = 0.4 nm. (b) Tip state: p z orbital, sample state: 
p z orbital, Rf, = 0.27 nm. (c) Tip state: d z 2 orbital, sample state: p z orbital: R c = 0.2 nm. (d) 
Tip state: f z s orbital, sample state: p z orbital, = 0.15 nm. 



p z - d z 2 : R = (5) 

1 + 4/ KZ 



p z — f z 3 : R = — ; — (6) 

y J 1 + 7/kz v ; 

where k is the decay constant of the tunneling current with a typical value of k — 1 A -1 . 

Fig. |3] shows four combinations of tip and sample atomic states. If tip and sample states 
are s-type, the radius of the atom image is simply the distance between tip- and sample 
atom. The higher the angular momentum, the smaller the radius of the atom images. 
Sharply peaked atom images are observed with tip states with a large angular momentum. 
It is evident from Fig. 01 that the angular orientation of tip and sample orbitals is important. 
Fig. 0] displays only tip orbitals oriented in ^-direction, while the two other p-orbital orien- 
tations p x , p y , the four other d-orbital orientations d^, d xz , d yz , d x 2_ y 2 and the six other 
f-orbital orientations can also contribute to the tunneling current. However, it can be shown 
that the matrix element between a p 2 -sample and p 2 , d z 2 and f z3 tip is much greater than 
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the matrix element between a p 2 -sample and the orbitals that are not aligned in z-direction. 

An experimental test of Eqs. (jUKl) can be performed when imaging a surface with known 
electronic states with a tip that is also characterized by a known state. On the Si (lll)-(7x7) 
surface, the sample states are well approximated by p z states. The tip states depend on the 
chemical identity of the tip atom. We have therefore looked for tip materials which exhibit p, 
d and f states. Si with an electronic configuration of [Ne]3s 2 2p 2 is a good choice for p-states, 
because the 3p states are the highest occupied orbitals in Si. Co with an electronic config- 
uration of [Ar]3d 7 4s 2 is expected to have d-symmetry and rare-earth elements are excellent 
candidates for f-type tip atoms. We have chosen Sm (electronic configuration [Xe]4f 6 6s 2 ) 
for an f-type tip, because it is one of the most stable rare-earth elements. For imaging with 
very small tip-sample distances, large forces are expected to occur and mechanical stiffness 
of the STM tips is important. Rare earth elements are interesting for STM studies, because 
in contrast to hydrogen and other lighter atoms, the 6s electrons are bonded more strongly 
to the nucleus than the 4f electrons. Thus, the 4f electrons participate in electronic conduc- 
tion and in the tunneling process from tip to sample. However, the average distance to the 
nucleus is larger for the 6s electrons than for the 4f electrons. Imaging Sm by STM implies 
that the outer 6s electron shell has to be penetrated - literally requiring to peek into the 
inner parts of an atom. 

Fig. El displays the series of dynamic STM images of p, d and f atoms, realizing the three 
cases presented in Fig. |U (b)-(d). In Fig. El (a), a Si front atom is imaged by the Si surface, 
in Fig. El (b) a Co atom and in Fig. El (c) a Sm atom. As revealed by this image, with 
increasing angular momentum of the tip states, the apex radii decrease indeed, very much 
as predicted in Eqs. (jHEJ). Figure El (c) is recorded with a tip made of pure Sm. The clear 
visibility of the rest atoms and the enhanced corrugation together with the very small apex 
radius of approximately 0.14 nm prove the angular confinement resulting from 4L,3 (/ = 3) 
tip states. Rest atoms have already been observed by force microscopy j^lSjQl- 

Fig. El (a) shows an experiment where a CoSm rare-earth magnet mounted on a quartz 
force sensor was used as a tunneling tip. The composition of the tip was CoeFeaSm, as de- 
termined by Energy Dispersive X-Ray (EDX) analysis. The magnetic field axis was oriented 
perpendicular to the sample surface. The eigenfrequency of the force sensor was /o = 19621 
Hz, and the positive frequency shift which was measured during imaging implies repulsive 
tip-sample forces during imaging. Qualitatively, the adatom images in Fig. El (a) can be 
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FIG. 5: (Color online) Dynamic STM images of the Si (lll)-(7x7) surface showing the influence 
of the angular momentum number I of the tip state to the apex radii. The images have been 
taken with one instrument, using similar imaging parameters but tips made of different materials. 
The spacing of the contour lines is 20 pm and their increasing density from (a) to (c) is caused 
by decreasing apex radii, (a) Image taken with a Si tip (a Fe tip suffered a mild collision with a 
Si sample, leading to a Fe tip covered with a Si cluster). Sample bias 300 mV, tunneling current 
lav = 300 pA with A = 0.1 nm. The tip atom is presumably in a p 2 -like sp 3 -state (1 = 1). The 
experimental apex radius is approximately 0.5 nm. (b) Image taken with a Co tip at a sample 
bias of 200 mV, tunneling current I av = 100 pA with A = 0.5 nm. The image is expected to 
originate from a convolution of a Co 3d 2 2 state (Z = 2) with the p z Si states. The apex radius is 
approximately 0.2 nm. Rest atoms are visible in the left half of the unit cell, (c) Image taken with 
a Sm tip (see also Fig. |2t). The rest atoms are clearly separated from the adatoms. The apex 
radius is only 0.14 nm. The tilt angle of the f z 3 state is presumably close to 0°. Other states (4f, 
6s) are possibly contributing to the image, too?" ^ 
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FIG. 6: (Color online) (a) Ultra-high resolution image of a Si (lll)-(7x7) surface, acquired by 

dynamic STM with a CoeFesSm tip oscillating with an amplitude of A = 0.5 nm at a sample bias 

voltage of -100 mV and an average tunneling current I av = 200 pA. The green arrow indicates an 

atomic defect (missing center-adatom) . (b) Schematic presentation of the current carrying atomic 

states leading to the observed image shown in a. The sp 3 silicon states are tunneling mainly into 

a Sm 4f^3 tip state tilted by a fixed angle, determined to be approximately 37°. (c) Experimental 

image of a single Si adatom imaged with a CoeFesSm tip. Average tunneling current I av = 1 nA, 

sample bias: -100 mV, amplitude A = 0.5 nm. (d) Calculated dynamic STM topography image 

for a silicon 3p 2 sample state and a Sm 4f z 3 tip state inclined 37° with respect to the z— axis. The 

calculation is based on the "Modified Bardeen Approach" j| in which the cantilever oscillation 

(amplitude A = 0.5 nm) is accounted for. Average tunneling current I av = 1 nA, sample bias: 

-100 mV. The color scales in (c) and (d) are identical, (e) Trace through maximum of image in 

x-direction (red line). The paraboloid fitted to the trace has an apex radius of R x = 0.12 nm. 

(f) Trace through maximum of image in y-direction (green line). The paraboloid fitted to the 

trace has an apex radius of R y = 0.15 nm. The black lines in (e) and (f) show the corresponding 
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cross-sections of the constant current surface at the main peak of the theoretical image. For clarity, 



explained with a p z -sample which tunnels into a f- tip state tilted by an angle G as shown 
in Fig. |U1 (b). Repulsive forces are expected to occur when the 6s shell in Sm has to be 
penetrated in order to allow for a large tunneling current. A comparison of this data with 
images taken with pure Sm, Co and Fe tips suggests that the tip atom in this image was Sm. 
The crescent surrounding the cusp is expected to originate from the upper hoop surrounding 
the club of a f z 3 orbital as shown in Fig. 13 (b). Several investigations were performed to 
check if these images are caused by experimental artifacts. Interference due to a multi-tip 
image has to be ruled out because of the presence of a single atomic surface defect indicated 
by a green arrow in Fig. |H1 (a) and because of the large depth of the corner hole image. 
Further, rotating the fast scanning direction and varying the scanning speed did not change 
the images, which rules out feedback artifacts. The cantilevers frequency shift and damping 
were recorded in parallel to the topography record. The frequency shift according to Fig. El 
(a) was positive, i.e. repulsive forces have been acting. The unusually sharp apex radii 
can not be explained by atomic relaxations due to tip-sample forces. While the radius of 
the atomic images decreases in STM when attractive tip-sample forces act 22j ]. imaging at 
repulsive forces leads to increasing radii of the atomic images, as shown experimentally in 
Figs. 15 and 16 in (l^ . In the experiments presented here, repulsive tip sample forces were 
acting which cause even an increase in the experimental apex radii. The damping signal 
did not show significant variations within the adatom image, which proves that a hysteretic 
lateral jump of the adatom or tip atom which possibly could have explained the sharp atomic 
peaks, did not occur [3]. In theory, fluctuations of the cantilevers amplitude could explain 
modulations in topography. However, the oscillation amplitude was kept constant by an 
amplitude control circuit and monitored, ruling out amplitude fluctuations as a source of 
the measured atom shapes. 

There is also quantitative theoretical evidence that the observed features are images of 
orbitals. Fig. |H1 (c) is a magnified view of an experimental single adatom image acquired with 
the CoSm tip in dynamic STM. Fig. El (d) shows a calculated image using the "Modified 
Bardeen Approach" which is based on quantum mechanical perturbation theory. For 
the many-electron-atoms silicon, cobalt, iron and samarium, " Slater- Type-Orbitals" 23^ 
were calculated as approximations for the electronic 3p and 4f states. The only free input 



parameter for the calculations is the electronic density of states at the Fermi level 
the tip, which was set to 0.4 electrons per atom and eV as a typical value for metals 
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FIG. 7: (Color online) Experimental dynamic STM image of the Si (lll)-(7x7) surface using a 
samarium tip. Sample bias -0.8 V, tunneling current I av = 50 pA with A = 0.8 nm. Two single 
atomic defects are present. 

The main result of these calculations is an excellent agreement between the experiment and 
theoretical results for Sm 4f 2 3 wavefunctions inclined 37° with respect to the axis normal 
to the silicon surface. Fig. |U] (e), (f) show contour lines in x and y direction for both 
experimental and calculated images. The calculation accounts for an oscillation of the 
cantilever with an amplitude of A = 0.5 nm in z— direction. In experiments as well as in 
calculations, the resulting image is insensitive to the amplitude A, as long as A is in the 
range of 0.3 - 1.5 nm. The difference between the experimental and theoretical images (a 
slight compression in z— direction of the experimental image compared to the theoretical 
image) is explained with atomic relaxations. Molecular dynamics calculations reveal that 
the repulsive interaction occuring at the lower turning point of the oscillation compresses 
the atom image in the z— direction, leading to an increased experimental atom radius. The 
excellent agreement between Fig. H (c) and Fig. |U] (d) implies, that in this experiment a pure 
4f 2 3 state has carried the major part of the tunneling current and contributions of other 
4f and possibly 6s states are negligible. Whether the high selectivity for the 4f 2 3 state is 
caused by the magnetic field generated by the permanent magnet or by crystal field splitting 
in the field of the neighbor atoms of the front atom remains to be studied. In standard STM 
experiments, much higher tunneling voltages than in our experiment are usually applied, 
probably resulting in a contribution of atomic states of more than one symmetry type (see 
Fig. El (a), (b). In this case, it is likely that the contributions of the various orbitals add to a 
roughly spherical symmetry of the image. As another important result of the calculations, 
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it was found that the oscillation of the tip favors a monotonic control signal for the z- 
control of the STM feedback. Therefore an oscillating tip may avoid a tip crash and enable 
lateral movement of the probe during the imaging process in cases, where imaging might 
be impossible using a static probe. For static tips, atomic resolution imaging in contact 
mode should be degraded because of instabilities of the tip atom. It was impossible to 
obtain atomically resolved images with a sample bias of only -100 mV when reducing the 
cantilever amplitude below 0.2 nm. Amplitudes considerably larger than the decay length 
of the tunneling current result in lower signal, because then electron tunneling can only 
occur for a small fraction of the oscillation cycle. The large spring constant of the cantilever 
used in our experiment assures stable operation at amplitudes in the range of interatomic 
distances. 

In Fig.[7[ dynamic STM data acquired with a samarium tip are shown. The special shapes 
of the atomic images, mainly with 4-fold symmetry, are supposed to be the consequence of a 
convolution of the sp 3 silicon states with different samarium states protruding from a single 
samarium atom. The special atomic shapes are repeated in the adjacent unit cells according 
to the periodicity of the Si(lll)-(7x7) surface. Here, an interaction with the rest atoms can 
not be detected. Single atomic defects in adjacent Si unit cells prove, that a single tip atom 
was responsible for the special atomic images. 

In summary, we have found a subatomically varying transition probability for the tun- 
neling process in STM, demonstrating the capability of dynamic STM to image structures 
within atoms caused by atomic orbitals. The observability of the substructures is attributed 
to the dynamic STM mode with a cantilever operated at amplitudes A in the range of 0.3 
nm < A < 1.5 nm and a small tunneling bias voltage. Combining dynamic STM with tip 
characterization tools like field ion microscopy 2^1 should further improve our understand- 
ing of the relation between atom shapes and tip states. Theoretical considerations link the 
experimental shape of the atoms to the atomic orbitals participating in the tunneling pro- 
cess. It is conceivable that chemical identification of the tip atom is possible with refined 
calculations of the tunneling current. 
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